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GENERALIZED AERODYNAMIC COEFFICIENT TABLE STORAGE, CHECK-OUT 
AND INTERPOLATION FOR AIRCRAFT SIMUALTION 

SUMMARY 

. For the last - three years the_ set of programs described in this paper has 
been used for rapidly introducing, checking out and very efficiently using 
aerodynamic tables in complex aircraft simulations. 

The preprocessor program reads in tables with different names and dimensions 
and stores them on disc storage according to the specified dimensions. The 
tables are read in from IBM cards in a format which is convenient to the engineer 
or person who reduces the data from the original graphs. During table processing, 
new auxiliary tables are generated which are required for table cataloging and 
for efficient interpolation. In addition, DIMENSION statements for the tables 
as well as READ statements are punched so that they may be used in other 
programs for readout of the data from disc without chance of programming 
errors. Besides punched card output, there is also a line printer output which 
consists of the properly labelled tables. For quick data checking graphical 

f 

output for all tables is provided in a separate program. 

Linear ^interpolation routines were written which required only onesearch per 
argument, independent of the number of lists for each argument. , The linear 
interpolation routines have been written in FORTRAN to interpolate 1 dimensional, 

2 dimensional, and 3 dimensional tables. 

Ihe programs described so far have been written to be used for aircraft simulation 

on the IBM 360. For piloted simulations of the EAI 8*1-00 computer the same 

< 2 . 

set of data is required. However, the data. input format is quite diffe^nt. . 

To save time and reduce operator errors a program was written that converts the 
format for the IBM 360 to a punched card deck in the EAI 8b-00 compatible format. 
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INTRODUCTION 


The programs that will he described in the following sections are designed to 

speed up the data storage, data check-out, and table look-up for aerodynamic 

/ /a 

coefficient tables used in aircraft automatic control simulations on the IBM 

360 and for piloted simulations on the EAI 8^00. 

On the 360 everytime a new aircraft or spacecraft had to be simulated new data 
processing subroutines had to be written, while the 8400 interpolation is not 
too efficient in its argument searches. The new table preprocessor program 
that will be described does not require any re-writing as function of the table 
names, table lengths, or table dimensions. The new interpolation program 

performs only one search per argument •_ 

independent of the number of argument lists required for the set of tables. 

This is further explained below. 

For an elaborate 6 degree of freedom aircraft simulation, up to 100 tables 
are required. These tables often have common arguments, for instance 'alpha' 
the angle of attack, but the list of argument values at which the function 
values are stored may be different depending on the shape of the functions which 
are approximated as linear segments. These lists are called break-point lists. 
During a simulation run each break-point list must be searched for the subscript 
of the break-point just below the actual argument value. Therefore, many 
break-point lists may have to be searched even though they concern the same 
argument. In our table look-up method this. disadvantage is over-come by generating 
master lists of all break-points and searching the master lists only. Hie 
subscript of the appropriate break-point for each argument of each table is 
. found in an auxiliary cross reference list. Hence, there is only one- search 
per argument. An inefficient search could also be avoided by expanding tables 
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by means of interpolation so that they all have the same break-point lists. 
If one assumes infinite storage, this is the most efficient method of table 
look-up. However, table expansion may increase the required storage by a 
factor of 4 to 8 so that a large amount of virtual memory and core swapping 
my be required on the IBM 360 which would make the program run inefficiently 
although the calculations would be reduced. In the programs that will be 
■described, the tables are kept at their original size at a cost of a small 
amount of computation. 


At present, data packages for aerodynamic tables for aircraft simulation are 
prepared separately for the IBM 360 and the EAI 8400. This requires two 

sets of data decks, each containing different types of data transcription 

■ \ 

and key punch operator errorsjand therefore, it requires two check-out periods 
Such check-out is very costly in computer and engineer's time. Therefore, a 
program was written that circumvents double data check-out by converting the 
data stored and checked out on the on the 360 to a deck of cards in the format 
that is required for the 8400 for aircraft simulation. 
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CONSIDERATIONS FOR THE PREPARATION OF TABLES 

Often aeronautical tables are generated from graphical data. The following 
remarks are intended to aid in speeding up the data reduction process and 
aid in obtaining an efficient table lookup for running simulation programs. 

The number of breakpoint lists should be kept small. This can be done by 
examining all tables with identical parameters and by selecting a set of 
parameter breakpoints that would permit faithful description of all curves 
with that set of breakpoints. Then select from this master list a reasonable 
number of subsets to describe individual curves with a minimum number of 
breakpoints. The maximum number of breakpoint lists .per argument is 20, 
but considerable storage and running time is saved for the simulation program 
if this number is kept smaller. Greatest care in reducing the number of break, 
point lists should be taken for the argument that appears most' often in the 
tables, usually the angle of attack. 

All entries of the tables must be filled. It often happens in graphical data 
that all curves on a graph do not span the full range of the parameter on the 
abscissa. It is thought that the engineer can estimate a better extrapolation 
for the purpose of filling the table than an automatic extrapolation program 
could. If no prior knowledge is available, or if the simualtion would be 
invalid if it entered the region of the graph for which no data exist, the 
engineer can insert extremely large function values, which would make the 
simulation fail in a very. obVious way. 

Effort is saved when the tables are constructed from the original graphs by 
using the same .convention as the automatic plotting program. Unless specified ' 
differently, for an individual table, the plotting program plots the horizontal 
breakpoint list as abscissa, the vertical breakpoint list as curve parameter, 
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and for a three dimensional table, it uses the third breakpoint list as fixed 

> / 

parameter for each individual graph, see figure 1. 

In naming the table it is useful to adhere' to a simple convention similar for 
all tables. A suggested -way is to add a T to the name of a function for the 

i 

tabulated values . Then, e.g. the table name for the lift coefficient CL 

would be CLT. 

The table parameters may -be named identical to the names given in the 
simulation program. The parameter names in the table processoriare used only 

as^phameric data. They aid in labelling the printouts and graphs and are 

'used in error messages in the simulation program. 
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Vert BPL 



iWU BPL HorU BPL 


(a) Graphical data. 



(b) Tabulated data. 


Relationship between breakpoint lists and standard plotting convention 






TABLE PRE-PROCESSOR 

READING IN THE TABLE DATA AND DESCRIPTION OF THE PROGRAM FLOW 
Tie table pre-processor consists of a main program and 5 subroutines all 
written in FORTRAN (for listings see appendix I). Tbeir operation will be 
explained by means of a small example. Figure 2 shows the data deck for the 
example, which includes short 1, 2, and 3 dimensional tables in arbitrary 

order. The cards are numbered and anotated for the purpose of discussion 

> 

only. 

The main program T READ, reads in the NAMELIST "Input", which contains only 
the variable IPUNCH. Table storage and card punching is controlled by IPUNCH 
when it is 'O' only line printer output occurs. It is advisable to run the 
program in this manner until one is reasonably sure of the correctness of the 
data format and table dimensions. Finally, set IPUNCH to '1' to store the 
data on disc and obtain a deck of DIMENSION and READ cards for all tables and 
auxiliary arrays. 

TREAD then calls subroutine TABCOM. TABCOM reads two key cards, (l) and 
(2), which contain the number of arguments for all tables to be read and the 
alpha-numerical names of all arguments . For reading each table TABCOM calls 
TABIN. Each table has a header card, (3)> (7)> ( 20 ). The header card contains, 
in order, a user's comment field, the table name, the names of the arguments, 
and the lengths of the individual breakpoint lists which are also the table 
dimensions. As can be seen, for the examples of the 2 and 3 dimensional 
tables, the tables are constructed in a conventional fashion with one breakpoint 
list written horizontally across the top, (9), (l5)j (2.1), and the other 
breakpoint, list, vertically on the left (first, column of numbers). .The ■ s 

breakpoint for the third dimension are given at the head of each 2-dimensional 
subtable, (8), (l4). This arrangement requires repetition of breakpoint lists 
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(2) ALPHA DA 


MACH 


no. of arguments 
argument names 


ALPHA 


Table I 
I di mensiana I 


iCL.DA 

■r«V.5 


ALPHA 


MACH 



0. 

5. 

10. 


-10. 

. '^rr - 

. 1 

. 1 


0. 

.3 

.2 

.1 


15. 

.501 

.302 

. 4 


,20 . 

.8 

.4 

.598 



0 . . 
.05 
.14 6 
.46 
.39 


5. 

.048 
. 1. 
.15 
.20 


10 . 
.05 
. 098 

.2 

.31 


Tatle 2 

3 d /me-nsio na / 



CLMiM 

DA 

ALPHA 

3 

; 3 

1 


-9. 

10 . 

18. 




0. 

1 0.00012. .2001 

.41. 




6 a 

.2 

• 6 

.3 




10. 

.4 

1.0 

1 .2 





Ta t /e 3 

2 J i mens/^Tia/ 


Figure 2. Card images with explanations of a sample set of input data to the 
pre-processor. (Blank cards were inserted for printout clarity only.) 


for 3-dimensional tables, but it ails the visual check of the data cards 

before processing, lie corresponding function values are located at the 

intersection of one element in the vertical breakpoint list and another element 
•in the horizontal breakpoint list. 


In addition to the above, TABCOM stores the function elements and saves the 
name of the table and the names of the arguments associated with the break- 
point sublists. Next, once per table, TABCOM calls TABOUT which is a Subroutine 
t writes out the table. The output for the example in Figure 2 is shown 
in Figure 3. TABCOM also compares the new breakpoint lists with those already 
stored and adds them to the already stored breakpoint lists, unless an identical 
breakpoint is already in storage. Furthermore, TABCOM assembles breakpoint 
identification lists which identify the required breakpoint lists as function 
Of the table number and first, second, and third table dimensions as appropriate’. 
After all tables have been read in, TABCOM prints out all breakpoint lists as 

shown in Figure 4(a), Control is then turned back to the main program. 


The last job to be performed by the table pre-processor is to generate the 
lists that are required for a more efficient table look-up. First the main 
program calls subroutine MLIST. MLIST generates a master list for each argument 
that occurs in the tables. (A master list is an ordered list of all values 
that occur in the sublists for one particular argument). These master lists 
are required for searching each argument that occurs in the tables. The master 
lists are printed out as shown in Figure 3(b). Finally, subroutine HBLIST is 
called. This subroutine generates cross reference lists which will allow table 
look-up by only searching the master list for each argument. After these 
lists are constructed they are printed out as shown in Figure 5 (a.) for the 

example, ,The NBPL lists define the corresponding element in' the breakpoint ^ 
list for each argument and each table. When the breakpoint list for a given 
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NUMBER OF VARIABLES 3 
4LPJ1A DA MACti 


/ 


table 

i 

CM 

( ALPHA ) 



l 

> ALPHA 

CM 



0.000000 

5.000000 

-0.200000 
O.fiOOQOQ 



10.000000 

0.800000 


TABLE 

2 

CLDA 

( ALPHA 

, DA , MACH ) 




MACH 

= 0.500000 


Al PH A 



0.000000 

5.000000 

10.000000 


DA 





-10.000000 





0.000000 

0.300000 

0.200000 

0.100000 


15.000000 

0.501000 

0.302000 

0.400000 


20.000000 

0.800000 

0.400000 ' 

0.598000 





- .-X 

- 

rV.V;V 

Figure 3. Printed tables from the 

preprocessor. 

. .. v' 1 ’ ‘ ’»*> ’ 


- 10 - 




MACH 


0.900000 


ALPHA 



0.000000 

5.000000 

10.000000 


DA 


-10.000000 

0.050000 

0.048000 

0.050000 

0.000000 

0. 146000 

0. 100000 

- 0.098000 

15.000000 

0.460000 

0.150000 

0.200000 

20.000000 

0.390000 

0.200000 

0.310000 






TABLE 

3 CL MN 

( DA 

, ALPHA ) 


DA 

-9.000000 

10.000000 

18.000000 

Ai PHA 

0.000000 

0.000120 

0.200100 

0.410000 

6.000000 

0.200000 

0.600000 

0.800000 

10 .000000 

0.400000 

1.000000 

1.200000 


(b) 

FIGURE 3 (concluded) . 
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ALPHA 

1 CM 0.00000 5.00000 10.00000 

3 CLMN 0.00 000 6.00000 10.00000 


DA 

2 CLDA * -10.00000 0.00000 

3 r.l M N -9.00000 10. 00000 


15.00000 20.00000 

18.00000 


MACH 

2 CLDA 0.50000 0.90000 


(a) Break-point lists, none repeated, with the table number and name where 

each list first occurred. 


MASTER LISTS 
4 ALPHA 

0.000 5.000 

6.000 

10.000 

7 DA 



-10.000 -9.000 

0.000 

10.000 15.000 18.000 20.000 

2 MACH 




0.500 0.900 


(b) Break-point master lists with number of elements and name of 

argument. / 

Figure 4. Auxiliary arrays from pre-processor. 



NBPL LISTS FOR 

ALPHA 



CM 

CLMN 

1 0.000 

1 

1 

2 5.000 

2 

1 

. .. 3 6.000 

2 

2 

4 10.000 

3 

3 


NBPL 

LISTS FOR 

DA 




CLDA 

CLMN 

1 

-10.000 

1 

-1 

2 

-9.000 

1 

1 

3 

0.000 

2 

1 

~-4 

10.000 

2 

~ 2 

5 

15.000 

3 

2 

6 

18.000 

3 

0 

7 

20.000 

4 

0 


NBPL LISTS FOR 

MACH 


CLDA 

1 0.500 

1 

2 0.900 

2 


(a) NBPl cross reference lists. 


CROSS REFERENCE LISTS 

J NTDIM NJ1 NJ 2 NJ3 NK1 NK2 NK3 NL1 NL 2 NL 3 

1 1 1 ; _1 3 

23 1 23 11 1 3 42 

3 2 2 1 2 2 3 3 


(b) 


■fcabfe 

Y»e. 


no. 

si i men. 


u — — / 

n,o. 

■f»r eac.K 


br«a|5-p®int li’s't no. leng'tK of 
for the, • ^Wen arj.j break-P®' lists 
1st, 3nd <jlme,n. 


Cross reference lists indicating the correct breakpoint list for each 
argument of each table. 


16 DIMENSION CARDS WERE PUNCHED. 


18 READ CARDS WERE PUNCHED. 

NOTE ON TSSt THESE CARD IMAGES ARE STORED IN THE DATA SET 

ASSIGNED TO FT07F001 BY A DDEF CARD. 

. SE E THE TABLE PROCESSOR USERS GUIDE. 

(c) Message giving number of access card images stored.. 

Figure 5. Final output from the pre-processor. 
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table does not cover the whole range of the master list, for points outside 
the range, a'-l' is inserted when the value is below the range and a 'O' 
is inserted in the NBPL list when the value is above the range. (See the 
DA NBPL list for table CLMN) . As will be shown later, this permits the inter- 
polation subroutines to determine, with a single IF statement; whether an 
argument is outside the range of the table. 

Other cross-reference lists that are 'required are shown In figure 5(t>). A't 'the 
beginning of the main program, under symbol explanation, the meaning of the 
important variables and tables are given. Sufficient comment cards are 
inserted into the programs to make them self-explanatory. Also given in the 
main program are the ranges of the subscripts, which are determined by the 
dimension statements presently in the program. The program as it stands 
can handle 25 different independent variables and a maximum of 20 different 
breakpoint lists per argument. It is of computational and storage advantage 
if the numbers of breakpoint lists per argument are kept to a minimum. The 
maximum length of a master list must not exceed 35 elements. These limits, of 
course, can be changed by simply changing the dimensions of the arrays. 

Since the pre-processor has stored all arrays efficiently, the remaining 
programs could not be written in complete generality. However, the only 
change that must be made as function of the type of tables to be processed 
is to take the dimension cards and read statements which are provided by the 
table pre-processor (see Figure 6) and insert them at the beginning of the 
main programs. Note that the READ cards must remain in the correct order. 



y 


nTMPWCTDM r M 

i ^ \ 

D-- 

_I 


DIMENSION CL DA 

( 3 , 4 1 2 ) 

D 

2 


DIMENSION CLMM 

( 3, 3) 

D 

3 


REAL-8 VARNMf 3) 


D 

4 


DIMENSION 8PL ( 3, 

2, 4) 

D 

5 


DIMENSION NLENTH ( 

3 r 2 ) 

D 

6 


DIMENSION N TABLE! 

3, 2) 

D 

7 


R E A L - 8 TNWE( 3) 


D 

8 


DIMENSION JTN ( 3) 


D 

9 


DIMENSION NTDIMf 

3) 

0 

10 


DIMENSION NBPL.f 3 

, 7, 2) 

D 

11 


DIMENSION N M AST' ( 

3) 

D 

12 


DIMENSION AM A ST ( 

7 i 3) 

D 

13 


DIMENSION NJ( 3, 

3) 

D 

14 


DIMENSION NK( 3, 

3 ) 

. D 

15 


DIMENSION NL( 3» 

3) 

D 

16 


READ ( 2 ) C M 



1 


READ (2) CLDA 



2 


READ (2) CLMM 



3 


READ (2) NUMTBL 



4 


READ (2) NOVARS 



5 


READ (2) VARNM 



6 


READ (2) R PL 



7 


READ (2) NLENTH 



8 


READ (2) MtABLE 



9 


READ (2) TNME 



10 


READ (2) JTN 



11 


READ (2) NT DIM 



12 


READ (2) MB PL 



13 


READ (2) NMAST 



14 


READ (2) AMAST 



15 


READ (2) NJ 



16 


READ (2) NK 



.17 


READ (2) NL 



18 

- 


Figure 6. Listing of the DIMENSION and READ cards provided by the pre-processor. 
These cards must be inserted into the graphics program, the 8400-card punch 
program, and the interpolation* program.. 
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• data, verification 


Although the tables are printed out in easily readable form, it is still 
tedious to compare the data with the original graphical aeronautical data. 

This process is simplified by graphical output of all tables that have been 
stored. A main plotting program and some subroutines are included in this 
table processing package to plot the tables stored by the pre-processor. 
Hardcopy 8x8- inch plots are produced on the S-C4020 cathode- ray- tube 
plotter . 

Sample plots of the example tables described previously are seen in figures 
7(a) through 7(g). Figure 7(a) shows a 1 dimensional table, CM. Figures 
7(b) and 7(c) show a 3 dimensional table, CLDA, which is plotted with the 
variables in order as stored by the pre-processor and has been automatically 
scaled on the stored data. Figures 7(d), 7(e) and 7(f) show the table CLDA 
plotted where the plotting program options are used. The plotting order has 
been changed so that the variable DA is now plotted as the abscissa. The 
3 variables can be plotted in any order the user desires. 

I 

Another option illustrated is specification of the minimum and maximum values 
of the function CLDA to be used on the graph. Figure 7(g) shows the format 
for a 2 dimensional table and further illustrates results of specifying the 
minimum and maximum graph values for the function CLMN and the variable DA. 
Instructions for specifying the options desired and the cards used for figure 7 
are described in Appendix VI. 

Listings of the programs to produce the data verification graphs are presented 
in Appendix II. Detailed instructions for their use are in Appendix VI. 
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•TABLE 2 CLDA 














Figure 7(dJ. Three dimensional table (part 1) showing a change in 
- plotting order from the graphs of figure 5(b) and (c). 
New minimum and maximum values of the function CLDA 



were also specified. 
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CHANGE FROM 36Q. JO. JWCH DATA CARDS 


As has been discussed in the introduction, a different set of punched table 
data cards is needed for the 8^00. In the 8^)0 the breakpoint lists must 
be entered separately and identified with different names for each breakpoint 
list although they concern the same argument. The function values for each 
table must be entered separately in a specific format which separates individual 
numbers by commas. Also, a slash is used as termination indicator for each 
table. The program that does the conversion from the 360 to the 8400 format 
consists of a main program (p8^00) and a subroutine (FTOB) which converts 
floating point numbers to a short alpha-numeric array representing the same 
number in the shortest possible format followed by a comma. 

As shown in Figure 8 for the example tables, the output of the program provides 
a listing and a deck of cards with the following functions. There is one 
THRUPUT card for each table, which gives the name of the table as well as the 
arguments in the order in which they are stored. «Each THRUPUT card printed 
image is followed by a line of comment which will aid in associating the names 
of the breakpoint lists on the THRUPUT cards with the proper breakpoint list. 

Secondly, all breakpoint lists are given and the name of the argument for each 
breakpoint list is defined. In addition, a comment card is included. Since 
each breakpoint list must have a different, name in the 8 E 00 program, a suggestion 
for accomplishing this is given at the head of the printout for this program, as 
shown in Figure 9* 

NONE OF THE RREAKPOINTLI STS PUNCHED IS DUPLICATED . 

ATTNTT I ALLY PUNCHED ALLlRE AKPOIMTL I STS ARE CALLED BY THE SAME NAME, 

■ SINCE WE CANNOT PREDICT WHAT THE NAMES OF THE INDIVIDUAL BREAKP ONTL I STS WILL BE. 

SUGGEST I0N = CHANGE CARDS AS PER EXAMPLE 

FR0M = - THRUPUT FAFB ( C J A TOA ) ' ' ~ 

TO = THRUPUT FAFB ( C JA I ,DA1 ) 

ACCORDING TO THE NUMBER OF 'T^_ BjREAKPQ_I NTL I ST FUR T HE GIVEN TABLE. • ; 

fTTTs NUMB E R IS GIVEN IN THE COMMENf LINE. 

THE VARBPT CARDS MUST BE CHANGED IN A SIMILAR MANNER 

FIGURE 9 . Printout from Data Conversion Program 


- 23 - 



THRUPUT CM (ALPHA 

COMMENT, BK PT LIST NO. 1 

THRUPUT CLDA (ALPHA 

) 

, DA 

, M ACH ) 

COMMENT, BK PT LIST NOS. I 
THRUPUT CLMN (DA 

COMMENT, BK PT LIST NOS. 2 

1 

, ALPHA 

2 

1 

) 

*THI S 
ALPHA 

IS BK PT 
VARBPT 

lst no. 1 for 

• 0 t 5 • ,10. 

VARIABLE 

ALPHA 

*THIS 

ALPHA 

IS BK PT 
VARBPT 

/ 

LST NO. 2 FOR 

VARIABLE 

ALPHA 

*THI S 

IS BK PT 

.0,6. ,10. 

/ 

LST NO. 1 FOR 

VARIABLE 

DA 

DA 

VARBPT 

-10. ,.0,15. ,20. 
/ 



*THI S 
DA 

IS BK PT 
VARBPT 

LST NO. 2 FOR 

—9 . ,10. ,18. 

VARIABLE 

DA 

*THI S 
MACH 

IS BK PT 
VARBPT 

/ 

lst NO. 1 for 

VARIABLE 

MACH 

CM 

POINTS 

.5, .9 
/ 



CLDA 

POINTS 

- . 2 , . 6 , . 8 
/ 



* 

MACH 

= 0.50000 

• 1,.1,»1 
• 3 , .2 , • 1 



'r 

MACH 

.501 ,.302, .4 
• 8 , • A- , • 59 8 
= 0.90000 



7D57T04TT, 705 : 

. 146, .1, .098 

• 46t • 1.5 , • 2 — — t — — 

CI.MN 

POINTS 

.39, .2 , .31 
/ 



.00012, .2001, .41 
.2,. 6, .8 
• 4,1. ,1.2 

/ v : 


Figure 8. Line printer output from the 8400 data conversion program. 
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One append s a number to the ; name of the' variable' for each breakpoint list' 
according to the number of the breakpoint list for the given variable. The 
comments automatically identify which breakpoint list number belongs to the 
variable for the given table. One must append this number to the name of 
the' variable in the THRUPUT card. We intentionally did not automatize this 
process, since the choice of the name for the breakpoint list is up to the 
user, and he may choose completely different names for the individual break- 
point lists. 

There are no changes required for card output of the function values them- 
selves. The proper comment cards are inserted for 3 dimensional tables. The 
function elements are ordered as specified in the THRUPUT statements with the 
fastest argument on the left and slower varying arguments successively to 
the right. The punched cards are identical to the line printer output, but 
in the 026 card punch format which is required for the EAI 8400 computer. 

Appendix VTI describes use of this program and shows a set of control cards 
for the TSS System. Again, the main program, shown in Appendix III must be 
recompiled for each card punch job because of -the DIMENSION and READ statements. 

With the availability of the data conversion program, the table pre-processor, 
and the plotting program it is hoped that in the future all function tables 
will be punched in the manner required for the pre-processor and coverted 
to the 8400 format by the conversion program. This should provide a faster 
data check-out, even if were not desired at present to put the aircraft 
simulation on the IBM 360. Also, hard copy graphical documentation of the 
tables that, are being used for . simulation are obtained. If, at a later date, . 

the program for an aircraft simulation is to be transferred to the 360, further 
benefit is derived from using the simple data card format for the pre-processor. 
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THE , INTERPOLATION SUBROUTINES 

i 

The interpolation subroutines are used to interpolate 1, 2 and 3 dimensional 
tables. Interpolation efficiency is increased by a search subroutine (SEARCH) 
which takes advantage of the cross-reference list NBPL, see Figure 5( a )* 

At the beginning of each programming cycle the master lists are searched by 
the subroutine SEARCH. This subroutine determines the subscript of the master 
list breakpoint if the nearest value is lower than the argument and calculates 
the proportionality constants PFAC for all breakpoint lists which are needed 
to interpolate the functions. 

[ Nearest loWt r~l 

PFAC - - Argument - L break- point J 

(Nearest kvgKarl^ |7?e<?r(»et /older”? 

re ak -/»£>• J jj> n t J 

First, however, the limits of the master lists are checked against the actual 
value of each argument. The argument is reset to the nearest limit when it 
is outside the limits and PFAC is set to zero. In this case a message is 
printed out indicating the facts. 


The acutal interpolation, of the tables is accomplished by three function 
subporgrams. FOUTID interpolates one dimensional tables; F0UT2D interpolates 
two dimensional table; and F0UT3D interpolates three dimensional tables. 

The three subprograms and interpolation equations are structurally similar. 


The one dimensional interpolation equation is: 

p - Ffct nearest: l<>k>c.rA + PFAC F /it nearest kMker) p/J't nearest Jew&rN I 
'*br« a *, t J v,trea kyv\ nt J \ #reaj--p e ; 0 t J J 

For the purpose of storage efficiency the subrpograms are not written in complete 

generality. The dimensions in the DIMENSION ■ statements must be changed to 

match, the array sizes required for a specific aircraft simulation. Since the 

* 

breakpoint lists of some tables do not cover the complete range of the master 
lists, each interpolation subprogram must check if the corresponding arguments 
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are within their limit. This is accomplished very quickly with the 

aid of the cross reference list NBPL, as explained in appendix VIII . 

Auguments outside the breakpoint list limit are reset to the appropriate limit for 
that table. 

For increased computational efficiency by a factor of four, the two 
and three dimensional table interpolation subprograms were written 
in one dimensional array lookup form. Listings of the search and 
interpolation subprograms are given in appendix IV. These sub- 
programs must be recompiled for each aircraft simulation to have the 
proper DIMENSION statements. Source programs may be obtained from 
the authors. Also, the READ cards from the pre-processor program have 
to be incorporated into the simulation to read in the tables. In 
appendix VIII an example is given of the programming required. 




CONCLUSIONS AND COMMENTS. 

A unified approach has been presented to table processing, checkout 
and table interpolation for aircraft simulation. Complete program 
listings are included in this report for documentation and to permit 
changes and additions by the user . The programs were specifically 
written with Ames computers and computer operation in mind, but portions 
of the programs could easily be adapted to other users' needs. 

Listings of the programs used in this table processjng package are given 
in Appendices I through IV. Detailed information for their use on 
the TSS system is given in appendices V through VIII. For further 
information on the control cards, users are referred to "A Guide 
to TSS/360 Batch Job Processing" and other TSS manuals available 
from the Computational Division library. > 

The computer programs comprising this table processing package have 
been compiled and stored in the authors' job library. This job lib- 
rary has been set up to be shared by all TSS users. Information 
on accessing it is given in Appendix V. 



APPENDIX I 


PRE-PROCESSOR PROGRAM LISTINGS 

This appendix contains listings of the computer programs used to 
process and store the aero tables. The programs (main and sub- 
routines) are in the following order: 

1- TREAD (main program) 

2. TABCOM 

3. TABIN 

4. TABOUT 

5. NBLIST 


6. MUST 
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APPENDIX II 

PLOTTING PROGRAM LISTINGS 


Data verification by plotting with the S-C 4020 plotter is done 
using the following programs: 


1 . 

Main 

2. 

PL0T1D 

3. 

PL0T23 

4. 

SCLIM 

5. 

FLIP2 

6. 

FLP213 

7. 

FLP321 

8. 

FLP132 

9. 

FLP231 

10. 

FLP312 


The subroutines PLOT! D and PL0T23 make use of subroutines in the 
Ames package of programs for the S-C 4020 plotter. 
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APPENDIX III 


EAI 8400 CONVERSION PROGRAM LISTINGS 


The following programs convert the tables to the proper format 
for input to the EAI 8400: 

1. Main 

2. FT0.B 
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APPENDIX IV 

INTERPOLATION PROGRAM LISTINGS 

The following programs perform table look-up and interpolation: 

1 . SEARCH 

2. F0UT1D 
, 3. F0UT2D 

4. F0UT3D 
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FUNCTION FOUTID(F,N) 

FOOT 10 INTERPOLATES A 1 DIMENSIONAL TABLE 
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FI = ‘F11 + PF2 *(F12-F11> 

F2 = F2 1+PF2 * ( F 2 2 - F 2 1 ) 

FQUT3D=F 1+PF3* ( F2-F1 ) 

RETURN 

END 



APPENDIX V 


USE OF THE PRE-PROCESSOR 


Genera] Information 

This appendix provides information to allow users to prepare their 
input cards and to process them with the table processing package 
described in this working paper. It is written to be a self-contained 
users guide. Consequently some of the information is similar to that 
given in the main body. 

The table processing package has been compiled and stored in a job 
library of the authors' TSS/360 library. A PERMIT command has been 
issued so that anyone can access the package on a read-only basis. 

Before the user can access the package, he must issue a SHARE command 
of the form 

SHARE DSNAME= MY .TABPROC ,USERID=FSMDNW1 ,OWNERDS=LIB. TABPROC 

where the underscored data set (job library) name is assigned by the 
user as desired. This SHARE command needs to be issued only once. It may 
be placed after the LOGON card for the first computer run. The user must 
include a DDEF card for his job library (e.g., MY. TABPROC) each time he 
makes a computer run using any part of the shared package. 

USE OF THE PRE-PROCESSOR 

The pre-processor set of routines is used to read in aerodynamic coefficient 
tables, process and store them, and generate auxiliary arrays for use in 
other parts of the table processing package. 



The first step in the process is to have the tables punched into 
IBM cards in the formats shown in figure V-l. This general format is 
similar to one commonly used when manually tabulating aero data from 
graphs. The card data fields are 8 columns wide, although such 
exceptions as table dimensions on the header cards should be noted. 
The table entries themselves are in floating point format and, for 
ease of card punching, are left justified. Note that everything is 
left justified in the fields except the integer values, which must 
be right justified. Floating point values must have a decimal point. 


The table header cards (cards (3), (7), (20) of figured (a) are 
set up as follows. The first field of 8 characters is a user's 
comment area. If desired for future reference, table numbers may be 
put here. For cataloginjpurposes , the Pre-processor assigns table 
numbers sequentially as the tables are read in. User's labeling of 
table numbers should follow this convention. The next fields contain 
the table name, up to 3 argument names as needed, and the break-point 
list lengths in the following order: 

Field Columns Contents 


9-16 

' Table Name 

17-24 

Horizontal break-point list name 

25-32 

Vertical break-point list name 

33-40 ' 

Third break-poin^ 1 ist name for 3rd tables 

41-44 

Length of horizontal break-point list 

45-48 

Length of vertical break-point list 

49-52 

Length of third break-point list 


Names are alphanumeric and are left justified. Lengths are integer 
and are rqight justified. Table names should be different from the 
function names they represent; for instance, CLT instead of CL. Break- 
point list names may be the same as the related variable names. 

































FIGURE V-lCb). Extended Format. 


If the 1st breakpoint list has more than 9 members this format applies. 
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The table entries are punched left justified into 8 column fields as 
in figured) (a). One dimensional tables are punched so that the first 
field contains the break-point (argument) value and the second field 
contains the corresponding function value, cards (4), (5), and (6). 

On 2 dimensional tables, the first card contains the horizontal break- 
point list starting in the second field, (2.1). Following cards have 
a vertical break-point value in the first field and corresponding 
function values across the card, (22), (23), (24). If a table has 
more than 9 horizontal break-point values and related function values, 
then overflow cards are made as in figured (b) . Three dimensional 
tables are made up similar to 2 dimensional ones except that each 
2 dimensional group is preceded by a card with the corresponding third 
break-point value in second field, cards (8) and (14) or figure'll (a). 

The format in figured must be strictly adhered to. Only a few 
checks are written into the program itself. Errors in column dis- 
placement are avoided for the floating point number entries by using 
a program control card (program 1) for the 029 card punch as shown 
in figured and by visually checking off-line listings of the data 
cards. For the table ID cards, one switches to program 2 of the pro- 
gram card. Program 2 is punched in such a manner that it is correct 
for 2 digit integers and for the alpha-numeric information for a three 
dimensional table. Hence, one must insert a leading blank or 0 for 
one digit integers and skip once or twice on tables with smaller 
dimensions. As a further restriction the horizontal break-point 
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FIGURE V- Z. 029 Program Control Card 



figure y-2(.a) 

lists must increase from left to right^(9 ), (15), (21); the vertical 
break-point lists must increase downward; and the break-point list 
elements for the third variable of a three dimensional table must in- 
crease in the order of appearance (8), (14). 




After the data cards are punched, one or more computer runs are 
made for initial checking of the tables. Figure Y-3 shows a sample 
\ deck of TSS/360 control cards for the computer run. 


CO 

(?-) 

(31 

(4) 

(5) 
<£) 


(7) 


LOGON US.E.RI.Q* JUS 1 22 JOHN DOF STOP 99 PH ,3999 

DDEF DDNAME = TABI. IB tDSQRG=VP , DSNAME = MY . TABPROC, OPT I ON= JOBL 1 R 

LDA D _TREAD$$ ~~~ __ 

'DDEF DDNAME=FT02F001 ,DSORG=VS ,DSNAME= TABLES 
DDEF DDNAME=FT07F001 » DSOR G=VS» DSNAME=DI MREAD 

CALL TREAD3>$ " 

£ INPUT 


i 


IPUNCH=1 , 
_LPJJ.NCH=_0_,_ 


SEND 


( 8 )( 


* 


L * 


TABLE DATA CARDS GO HERE 

,</ +>,- V*. J. sir >1. vt. O. X X J. J, 

'(* *i' v v *r v r v '<■* '•* v t ¥ v r 


(?) j 
(/o ; 


SEND 

LOGOFF 


J 


Figure V*3. TSS/360 control cards for computer 
run using the Pre-processor 

Unerlined words should be changed to conform to the user's naming 
conventions. Card (1) is the standard users LOGON, card. Card 

J 

(2) defines the user's job library which contains the table processing 
package. The data set name (DSNAME) should be identical to that 
specified by the user in the SHARE command. The data set name in 


-? 0 - 



card (4) specifies where the processed tables and related data are 
to be stored. Card (5) names the data set where the DIMENSION and 
READ card images a^red, The second of the NAMELIST "INPUT" 
cards (7) determines whether the tables and the DIMENSION and 
READ card images are to be stored (IPUNCH-1) or not stored 
(IPUNCH=0). The user's table data cards go at (8). The SEND 
at card (9) must be present to signal the end of the data cards. 


It is suggested that off-line listings of the data cards be examined 
for such things as blank cards, duplicate or missing cards, and the 
use of wrong columns. Further, it is suggested that the user make 
computer runs with IPUNCiM) to conserve computer time until all of 
his tables and the auxiliary arrays are printed out correctly. Some 
errors which may cause premature run termination include missing 
cards in a table and incorrect dimensions on a table header card -____ 
After all tables have been printed out, they should be examined 
careful ly for errors. Error message printouts should be noted and 
their causes remedied. Causes of errors can be determined by examinati 

of the printout and the data cards. 

After all apparent errors are corrected, another computer run should 
be made with IPUNCH=1 to store the tables and the DIMENSION and 
READ card images in their respective data sets. The DIMENSION and 
READ cards can be punched using the,\Command 

PUNCH DSNAME= DIMREAD ,STARTN0=1 ,ENDN0-80 
where the data set name is the same as specified on the DDEF card (5) 


of figure V-3. A printed listing of the DIMENSION and READ statements may 
be obtained by using the TSS command 

PRINT DSNAME = DIMREAD, PRTSP = 1 , 

STATION = RMT02, 

/ 

where RMT02 is for RJE station 2. These commands can be put before the 
LOGOFF card (10), run as a separate batch job, or entered conversationally 
from a terminal. The DIMENSION and READ cards should have interpreted and 
separated according to type for use in the plotting main program. They are 
numbered to aid in keeping the READ cards in order. 
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APPENDIX VI 

USE OF THE PLOTTING PROGRAM 

The plotting program is written to provide hardcopy plots of 1,2, or 
3 dimensional -tables on the S-C 4020 plotter. All necessary sub- 
routines are shared in the TSS job library described in Appendix V. 

The main program, however, is data dependent and must be changed by 
the user to include DIMENSION and READ cards for his stored tables. 

A copy of the main program may be obtained from the anthors. The 
required DIMENSION and READ cards are obtained during use of the 
pre-processor as described in Appendix V. Comment, cards in the main 
program indicate where these cards are to be placed. The main program 
may be compiled separately or during the plot job. 

The plotting program provides the user with several options for 
plotting the tables. The options and data cards to specify them are 
described with reference to figure VI-1. Card (1) transmits information 
to identify the plot job and to specify the user's tape number. Columns 
1 to 6&~are printed on the first and last frame of the plot job. 

Columns 69 to .72 are printed on the upper right hand corner of each 
plot frame. The user must have a 7 track magnetic tape assigned to 
him for this purpose. The plot option cards, (2) through (6), specify 
the table to be plotted, whether the variable order for 2 or 3 dimensional 
tables to be plotted differs from the order in storage and, if so, what 
order is desired. The other option specifies special scaling limits desired 
Special plot instructions can be largely eliminated by following certain 

corrections for the table construction as explained in the section 
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Card 
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Users 
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r\aw,£^ w\a / / S'ta } ^ of 
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Uob 

It 
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Kse.1 V=. 
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_ @@ (g 

@ -Us 

CO 

JOHN DOE 

STOP 99 PH 3999 PLOT 

JOB T 999 

X99 

99999 


(/£?) Plot f fl/'e'w+ ."fl cat / *k\. Card (must te ^• rs *’ data czrd). 


Plof Scaling OptUv\s 



(i) T.l <e p/o~i t?pt\cx\ cards. 


figTjr^ vi-1. example of user supplied information for table plotting 
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"Considerations for the Preparation of Tables". 


A logical variable called IORDER can be set to either -1, 0, or +1 
(column 8). If IORDER is equal to -1, all of the tables are plotted 
in the order in which they have been stored on the disc by the table 
pre-processor, see figure VI-2. In this case variable 1 is on the 
abscissa, variable 2 is always plotted as the fixed parameter for a 
given curve, and variable 3 is always plotted as fixed parameter for 
each set of curves for 3-dimensional tables. The scaling is automatic. 

If only a few tables are to be plotted from all the tables that are 
stored, IORDER may be set equa] to 0. In this case more specifications 
must be given. The second number (right justified column 16) will be 
the table number. As an optional feature, the next two numbers define 
the maximum and/or minimum value of the function to be plotted on the 
ordinate (columns 41 and 49) and the last two values designate the 

minimum and maximum of the first argument that is plotted as the 

abscissa (columns 57 and 65). Omission of these values results in 
automatic scaling of the data points as stored... No scaling, of course, 
is needed on the second and third arguments. As mentioned before, the 
arguments may not have been stored in the order as one wants to plot 
them out. In this case the logical variable IORDER is set equal to 1. 

In addition to the information that is needed when IORDER equals 0, 
one also must read in the order of the variables to be plotted by giving 
the names of the 2 or 3 variables in the proper order (columns 17,25, and 
33 .as._oji-xard._M) ) . 
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FIGURE VI-2. RELATIONSHIPS BETWEEN BREAKPOINT LISTS AND STANDARD 


PLOTTING CONVENTION 








The TSS control cards required for the plotting computer run are shown 
in figure VI-3. 




(2K 


• LOGON 'user ID T JOB r £9 (7) JOHN DOE ST OP 3L9 PH .3999. 
ERASE SOURCE. MYP LOT S £ ? FTN MYPLOT SSL* I SD = N 

»■ - ' mmfzwmf / 


Supplied main program w i th user ' s dimension and read cards inserted 


GOES HERE FOR COMPILATION 


■ *** ^ "Jr Vr . ^r ” I 

V, *|' 'l' *>' «yx '•,-x 'I' -v 'V *T“ 'J' 'f K' n' 'T 'T* 't* 'T 'I* ' 


. V. .1. S>/ xA. U< »*> 


(i) 

Of) 

(I) 

((>) 

(7) j 

( 8 ) 

ml' 


DDEF ddname= tablib ,dsorg=vp »dsname =my.tabprog » opti on = joblib 

AMES SC4020 

load myplotss 

DDEF DDNAME = FT02F001, DSORG=VS,DSNAME = TABLJES t 

MTMSG 

PLEASE OBTAIN" 9999 9 '7-T RAC K WITH RING 
CALL MYPLOT $ S 

_JL_ (J, -JL. .1 _J_ .t. U. xf. J. xV -J-. xA. xl. .V X*. xA. x V xA» V** "4 xC >V X*. xl> 'U X*. 

5^ 5,C ^Jc ^x w *r *r *I* 'i' *>* v 'I' ■i* *|* *»* "i" *y* v "p w *r 'i •i' *i* *r 


( 10 ) 

00 


* PLOT OPTION CARDS GO HERE 

J/ -£? -*- 

-r v t 

'jl li. xCxL xO x 1 / J. xC xC x>. x r . x*/ .1. xC x*/xC A A A A A A -A A A A A A A 
2,; J,C 5,C -r- V -v* -V -i- -v -i- -r '»• *v -r W v*r *i*¥ 'i*v *»* 

?END 

LOGOFF 

Figure VI-3. Control Cards for Plotting Job. 

Underlined parts of the cards are to be changed by the user as necessary. 
An estimate of time required, for the LOGON card, can make based on the 
authors’ experience in plotting the number- of tables below during 
1 job run. 

No. of Tables Time (sec) 


3 

32 


15 

132 


Use of the plotting options affects the time used, so it is suggested 
that 10 - 30% additional time be specified to avoid premature termination 
of the job. Portions of a job cannot be saved if the tape writhing job is 
terminated due to underestimation of time. The plotting main program, 
with the user's DIMENSION and READ cards from the pre-processor, are 
put at (2) with the FTN card for .compilation. The DDEF card (3) specifies 
the user's job library which contains the shared program modules. Card (4)' 
provides access to the Ames library of plotting routines. The LOAD card (5) 


must specify the main program named at (2). The user's stored tables 
to be plotted are specified by card (6). Card (7) specifies the user's 
7-track magnetic tape. Card (8) starts writing the plotting tape. 

The plot identification card and plot option cards described previously 
go at (9). Card (10) signals the end of the data cards. The printer 
listing received after the job run will specify which tables were 
plotted and the number of plot frames made. 

The plot run described above generates plotting data on the user's 
magnetic tape. To get hardcopy plots made, the user must send a card, 
as in figure VI-4, to the operations room, N233. The card is self- 
explanatory. Off-line production of the plots from magnetic tape 
takes approximately 24 hours. 


SC-402O PLOTTER 
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I 

! Tttv /rw 

i 

, 
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MICRO FRAME OFF 


OFF 


ON 


OFF 


Figure VI-4. Request for SC 4020 Plots. 



APPENDIX VII 

USE OF THE EAI 8400 DATA CONVERSION PROGRAM 


i 

This appendix describes the use of the programs which convert the tables 

stored on disc by the Pre-processor to a punched card format which 

can be read directly by the EAI 8400 computer. Only a small amount 
* 

of manual work, as described, is necessary on the punched cards. 

The main program is data dependent and must be recompiled by the user. 

A card deck copy of it may be obtained from the authors. A set of the 
DIMENSION and READ cards produced by the Pre-processor is inserted 
in the main program at the place denoted by comment cards. The main 
program may then be compiled alone or as part of the computer run 
to punch "8400" cards. 



(¥> 

(5) 

00 

( 7 ) 



LOGON USERID * JOB >99' "JOHN DOE " STOP 99. PH 3999 

ERASE SOURCE , MY P84$$ FTN MYP84SS* ISD=N 


SUPPLIED MAIN PROGRAM WITH USER'S DIMENSION AND READ CARDS INSERTED 
GOES HERE FOR COMPILATION 


. J. si. O, J. , 


DDEF DDNAME = TABL I B , DSORG =VP * DSNAME = MY , T ABP R PC , OPTI ON= JOBL I B 
LOAD M YP 8.4.$ $ 

DDEF DDNAME=FT02F001 ,DSORG= VS, DSN AME= TABLES 
DDEF DDN AMF = F T0 7 F QQ 1 «PS0RG=VS,0SNAME =CRD84 

C A LL MYP84S$ ~ • 

PUNCH DSNAME= CRD84, STARTNQ=1*ENDN0=80 
LOGOFF 


Figure VII - 1 . Control cards for punching data cards for EAI 8400. 



The TSS control cards used compile the main program and to execute it 
are shown in Figure VII-1. Underlined parts are to be changed by the 
user as necessary. Card (1) is the user's LOGON card. Card (2) 
is for compilation of the main program which is placed at (3). Its 
module name may be chosen by the user and used on cards (2), (5), and 
(8). The DDEF card (4) specifies the user's job library which contains 
the shared program modules. Card (6) specifies the user's tables 
which were stored by the Pre-processor. Card (7) specifies the 

dataset where card images of the converted tables are^to be stored. These 

cards are punched by card (9), which may be used here, run separately, 
or entered from the conversational terminal. 

The printed output from the computer run gives card image printouts 
which should be examined by the user. All breakpoint and table 
values punched out will be between 0.0001 and 9999.9. If a table entry 
should be outside this range, the number will be replaced by a series 
of slashes (/////.), which can be easily recognized and can be replaced 
by the proper number. Further information about changing the punched 
cards to obtain a running deck for the 8400 program is given in the 
main body of the report. 


-too- 



APPENDIX VIII 


THE APPLICATION OF THE TABLE LOOKUP PROGRAMS 

The use of the table lookup programs is illustrated by means of an 
example. The tables to be interpolated are those saved and plotted 
previously. The sample program is shown in figure VI 1 1 - 1 . First, the 
DIMENSION and READ cards from the table pre-processor must be inserted 
in the simulation program to make the tables and auxiliary arrays 
available to the simulation program (cards 500 to 3800, not shown 
on the figure). In addition, the common blocks "ARGS" and "CATLOG" 
must be put into the user's simulation program. Before calling 
the SEARCH routine, the table arguments must be placed in the proper 
order into the array ARG (cards 4400 to 4600). SEARCH resets the value 
of any member of the array ARG to the closest tabulated limit if 
the argument was outside the tabulated range of the master list. A 
statement is then printed to this effect. Figure V I I I - 2 lines 9 and 
10 are examples of this output. The lower limit for APLHA in the tables 
was 0. and the upper limit for MACH was 0.9. Note that SEARCH does 
not reset the actueal arguments ALPHA, MACH and DA as shown in line 
12 of figure VII 1-2. Other computations in the simulation program 
would therefore be based on the actual values of the arguments. This 

non-reset of actual arguments feature of the simulation can be traded 

/ 

for a small amount of increased computational efficiency by putting 
the actual argument list in the labelled common block ARGS in the proper 
order: COMMON/ARGS/ ALPHA, DA, MACH 




0 . i 


R F A L MACH . rfV r ~ — tt 

CP- MOM /ARCS/ ARP(^) r r T r-TqTi\T^TrS7T^iX^rTP¥'^ ^ ’ ' W " 

— r ;n^rn^---7i:TrrrnnTTFPfc program . .■ : 

„ T ;.^unN REAO CAROS ^^^LCtSb __ ...... 

X * - ■ I : 


2 00 


u; m'i 

;•! 'Sl^Io'tSa, ^.SO O) alpha, o^h ach 

ARR4¥ IM THE 56 


M„p ORDER A S FPTFREn I^TP 


JHE DATA S i OR AG 6 PR OGRAh 
-7TRCi”1r) — ^-GVtTi^A 
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/>. R G ( 3 ) = MACH 


HTST — SnRS'CR'TpT'S 


42 DO 
4300 C 
— ip-iifflfj 
4500 

4 40 0 - T - 

-^• 7 -oo---tr--F-tpn-^^^P^ r 

' 4R00 CALL SEARCH 

l 4900 C 
-5TTTTCT CT^ 

5100 
! 5200 
^3W 
• 5400 
5500 

5 700 1 'OCHA = » F 6 • 3 » / U-L 

5 800 CO T O 100 

~59 0 0 SOT) — S'TCrP - 

6000 ■ R,v '° 


— AND — CHiTFITTF FTCCP DR ' 


T •TM 7 rrD'-rY'''F ACTCP S 


TE THE TABLES 
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FCTUTZOTCr^ 73'1 


7TPTr"FI ' !•' H L K 
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c F R LIE T c T 1 1” 1 ? I R A L P H a , D & ■ j CH , C w ft . F1 L£Z- tt — -r ran 
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Figure VIII-1. Sample interpolation program 


ALPHA = 5.100 


da = 0.100 


G A C H 


0.510 


CFT'r~= o .7ZFT 

CLO = 0_._1.96 

POR = 0.347 


ZTCPHA 


nr=- — 7'. 05 0 nir^ 


ZTTiTO'n maUFT' - o . 7 0 0 




C 0 A - 0.120 



The function values are obtained when needed by function calls as shown; 
on cards 5100 to 5300 in figure VIII-1. /Lines! to 8 on the next 
figure show the argument values and results from the table lookup 
when all arguments are within the proper limits. If an argument is 
outside the proper limits of a specific table and within the range of 
the master list, a call to an FOUTxD function will result in an output 
as shown on line !!,, figure V 1 1 1 - 2 and the function value will be set 
to the proper sublist limit as indicated on line 14. 

'T"V ■*“ ‘ 



The structure of the NBPL list makes the above check very efficient. 

The appropriate NBPL element is tested with an algebraic IF statement. 
When the NBPL element, which normally determines the subscript for the 
breakpoint list, is negative, the argument is below the range of the 
breakpoint list for the given table. When the NBPL is 'O' the argument 
is above the range of the breakpoint list. In either case the follow- 
ing information is printed out: the number and nam^e of the table 

and the name ond the value of the argument which exceeds the limits 
of the breakpoint list. 

In some cases, it may be desirable to suppress the above error messages. 
To do this one can add a two value dummy table, which will not be used 
in the. actual program, but which serves to establish a breakpoint 
list for the parameter in question that spans the whole range of the 
parameter. An example might be the parameter 'altitude' for ground 
effect calculations. A second method to remedy the situation which 
is not quite as efficient, is to use logic statements in the running 


-10 2 - ; 


program to permit ARC(J) for the given parameter to take on values 
only within the specified limits. 

The table lookup routines are data dependent and must therefore 
be recompiled by the user after insertion of the proper DIMENSION 
statements. Card deck copies of these subprograms-SEARCH, FOUTID, 
F0UT2D, F0UT3D-may be obtained from the authors. 


All 4 of these subprograms have identical COMMON blocks and DIMENSION 
statements, as shown in figure V 1 1 1 - 3 . Most of the DIMENSION statements 
are provided by the pre-processor as described in appendix V. The 
sizes of the arrays ARG, LSUB, and PFAC are dependent on the sizes 
of VARNM and BPL as shown by the arrows in figure VI 1 1-3. The 
values of the array LSUB must be initialized to 1 at load-time. 

This may be done in a BLOCK DATA subprogram as in figure VI 1 1-4 
or by a DO loop in the user's initialization program. 
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Figure VII 1-3. COMMON and DIMENSION statements used in table lookup 
subprograms SEARCH, F0UT1D, F0UT2D, F0UT3D. 
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Figure VI 1 1 -4. Initialization of LSUB array. 
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